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This  discussion  compares  incident  and  transferred  spectral  results  for  gages  located 
at  Burns  Harbor,  Indiana.  Comparisons  are  made  by  examining  the  differences  in  the 
incident  and  transferred  harbor  energy  spectrums.  Spectral  analysis  allows  the  energy 
of  the  total  wave  record  to  be  broken  down  into  discrete  frequency  bands.  Energy  inside 
and  outside  the  harbors  may  then  be  compared  and  a  transfer  factor  for  each  discrete 
frequency  can  be  determined. 

Incident  and  harbor  wave  data  were  collected  within  and  offshore  of  the  breakwater 
at  Burns  Harbor,  Indiana.  The  purpose  of  this  data  collection  effort  was  to  determine 
characteristics  of  the  rubble  mound  breakwater  located  there.  Wave  records  were  col¬ 
lected  hourly  using  subsurface  pressure  sensors.  The  sample  rate  for  these  sensors  was 
1  Hz  and  the  burst  length  was  2048  seconds. 

The  analysis  utilized  the  Welch,  [  I  ],  spectral  analysis  method  with  50%  overlapping 
segments.  Since  the  raw  time  series  were  obtained  using  sub-surface  systems,  a  depth 
determined  high  frequency  cutoff  was  applied.  The  averaged  co-and  quad-spectra  from 
each  analyzed  record  were  used  to  calculate  significant  wave  height  (Hmo),  peak  period 
( Tp ),  and  mean  wave  direction  at  Tp  (Dp)  and  energy  spectrums. 

The  time  period  6  January  2003  0700  thru  7  January  2003  0  was  selected  for  de¬ 
tailed  analyses.  Figure  1,  top  plot,  shows  six  successive  hourly  energy  spectrums  for 
IN001,  located  in  the  open  lake.  The  lake  is  relatively  flat  at  0700  with  Hm0  0.53  m 
and  Tp  about  4.3  seconds.  Over  the  next  5  hours  the  Hmo  builds  to  1.86  m  and  the 
Tp  moves  to  6.6  seconds.  The  energy  at  the  long  period  end  of  the  spectrum  is  very 
small.  Examining  simultaneous  analysis  results  from  within  the  harbor  provides  a  de¬ 
scription  of  the  breakwater’s  performance.  Figure  1,  bottom  plot,  shows  spectrums  for 
the  same  six  hours  at  the  inside  gage,  IN002.  There  is  almost  no  energy  at  0700  GMT 
1/6/03,  with  Hmo  0.05m  and  Tp  4.5  sec.  In  subsequent  hours,  bimodal  spectrums  de¬ 
velop  as  more  and  more  energy  is  transferred  through  the  breakwater,  with  peaks  >20 
seconds  and  between  6  and  10  seconds.  The  magnitudes  of  these  peaks  are  small  when 
compared  to  energy  for  those  frequencies  outside  the  harbor. 

To  provide  a  more  direct  comparison  of  incident  and  transferred  energy,  a  transfer 
coefficient  (xfer)  can  be  calculated  by  dividing  the  transferred  energy  at  each  frequency 
by  the  corresponding  outside  energy,  eqn.  1. 
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where  E,-,  f  is  the  energy  per  frequency  transferred  inside  and  Eaf  is  the  incident 
energy  per  frequency  from  outside  the  breakwater. 

If  the  breakwater  is  working  well,  the  overall  transfer  rate  should  be  low.  Figure  4, 
bottom  plot,  shows  that  for  the  initial  hours  of  the  storm,  the  xfer  is  less  than  .15  for  all 
wave  periods  <10  seconds  for  IN002.  This  means  less  than  15%  of  the  outside  energy 
is  making  it  into  the  harbor  at  those  frequencies.  The  long  period  xfer  is  <1.0  which 
indicates  that  the  outside  long  period  energy  (>20  seconds)  is  still  more  than  that  on 
the  inside. 

After  the  first  six  hours  of  the  storm  the  Hm o  stays  around  2  meters  for  the  next  8 
hours  and  the  Tp  has  shifted  from  6  to  8  seconds  at  IN001  (Figures  2&3,  top  plots).  Lit¬ 
tle  incident  energy  is  evident  for  wave  periods  >20  seconds.  On  the  inside  at  IN002, 
figures  2&3,bottom  plots,  show  that  the  bimodal  spectrums  have  also  grown.  Since 
long  period  harbor  oscillations  cause  ship  movement,  the  energy  for  periods  >20  sec¬ 
onds  is  of  particular  interest.  There  appears  to  be  a  relationship  between  the  total 
energy  of  the  spectrum  and  energy  >20  second. 

Transfer  factors  for  1/6/03  1300  to  1/7/03  000  GMT  are  shown  in  figures  5&6, 
bottom  plots,  for  IN002.  For  some  of  the  periods  >  20  seconds,  the  transfer  factor  is 
greater  than  1 .0.  Another  observation  from  these  plots  is  that  the  transfer  factors  for 
waves  <10  seconds  appear  to  be  of  the  same  magnitude  throughout  the  storm. 

Figure  7  shows  the  directional  wave  statistics  for  January  2003.  This  plot  is  an 
example  of  the  wave  climate  for  IN001  and  shows  that  there  were  numerous  times 
during  the  month  that  the  Hm0  exceeded  0.5  meter.  In  almost  all  of  these  cases,  the  Tp 
was  greater  than  5  seconds.  Figure  8,  bottom  plot,  shows  the  mean  transfer  rates  for 
each  frequency  at  IN002  for  outside  wave  records  with  Tp  >5  seconds.  The  dashed 
line  represents  a  one  standard  deviation  band.  On  average,  the  breakwater  reduces  the 
energy  transferred  for  all  frequencies.  As  expected,  the  breakwater  performs  best  for 
short  period  waves.  Figure  9  shows  the  short  period  portion  of  the  Figure  8.  Mean 
transfer  rates  for  waves  between  5  and  10  seconds  is  less  than  5%. 

Portions  of  the  breakwater  has  undergone  repair.  Figures  8&9,top  plots,  illustrate 
mean  transfer  rates  for  a  gage  placed  in  front  of  a  repaired  section.  Energy  is,  on  aver¬ 
age,  reduced  for  all  frequencies.  In  particular,  energy  below  10  seconds  is  reduced  by 
50%  when  compared  to  results  from  a  gage  placed  in  front  of  an  non-repaired  section 
(Figures  8&9,  bottom  plots). 

For  more  information,  contact:  James  P.  McKinney  or  Margaret  A.  Sabol,  CEERD- 
HC-SO 
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Figure  1:  Top  Plot:  Six  successive  hourly  energy  spectrums  for  the  incident  gage, 
IN001,  located  in  the  open  lake  and  the  harbor  gage,  IN002,  bottom  plot. 


3 


I  NOD  1 


9.00 

7.20 

5.40 

3.60 

1.80 

0.00 


0.0 


0.1 


0.2 


igSEQUENCY ,  HZ 


0,3 


Hmo,  M  T  ,  Sffi  DEPTH.,  M 


o  = 

6-JRN-20Q3 

13  s  0 

0 

0 

1.90 

CJ) 

'vl 

12.8 

A  = 

6-JRN-20Q3 

1 4  s  0 

0 

0 

1.96 

7.1 

12.9 

+  = 

6-JRN-20Q3 

15  s  0 

0 

0 

2.09 

■"-J 

CO 

12.8 

X  = 

6-JRN-20Q3 

16  s  0 

0 

0 

2.20 

t-o 

CO 

12.8 

o  = 

6-JRN-20Q3 

17:0 

0 

0 

2.14 

■"-J 

CO 

12.8 

V  = 

6-JRN-20Q3 

18:0 

0 

0 

1.99 

CD 

CO 

12.8 

FREQUENCY,  HZ 


20 . 0 

10.0 

6 . 57 

5.00 

4.00 

PERIOD,  SEC 

Hmo>  M 

T  ,  SEC 

DEPTH,  M 

O  = 

6-JRN-2003 

13 

0 

0 

0 

0.18 

7.1 

7.8 

A  = 

6-JRN-2003 

14 

0 

0 

0 

0.20 

7.3 

7.9 

+  = 

6-JRN-2003 

15 

0 

0 

0 

0.22 

7.8 

7.8 

X  = 

6-JRN-2003 

16 

0 

0 

0 

0.24 

7.5 

7.8 

O  = 

6-JRN-2003 

17 

0 

0 

0 

0.24 

8.5 

7.8 

V  = 

6-JRN-2003 

18 

0 

0 

0 

0.23 

8.5 

7.8 

Figure  2:  Top  Plot:  The  next  six  successive  hourly  energy  spectrums  for  the  incident 
gage,  IN001,  located  in  the  open  lake.  Bottom  Plot:  Corresponding  inside  spectrum 
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Figure  3:  Top  Plot:  The  next  six  successive  hourly  energy  spectrums  for  the  incident 
gage,  IN001,  located  in  the  open  lake.  Bottom  Plot:  Corresponding  inside  spectrums 
for  IN002. 
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Figure  4:  Transfer  rates  by  frequency  for  the  first  6  hours  of  the  storm 
and  IN003. 
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Figure  6:  The  next  six  successive  hours  of  energy  transfer  values  by 
IN002  and  IN003. 
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Figure  7:  The  incident  directional  wave  statistics  for  January  2003 
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Figure  8:  Bottom  Plot:  The  mean  transfer  rates  for  each  frequency  at  IN002,  for 
incident  wave  records  with  Tp>  5  seconds.  The  dashed  line  represents  a  one  standard 
deviation  band.  Top  Plot:  The  mean  transfer  rates  for  IN003. 
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Figure  9:  Bottom  Plot:  The  mean  transfer  rates  for  each  frequency,  <10  seconds,  at 
IN002,  for  incident  wave  records  with  Tp  >5  seconds.  The  dashed  line  represents  a 
one  standard  deviation  band.  Top  Plot:  IN003 


